www.gusucode.com > matlab从零到进阶程序与数据 > matlab从零到进阶程序与数据/第18章 参数估计与假设检验/Chapter18.m
%-------------------------------------------------------------------------- % 第18章 参数估计与假设检验 %-------------------------------------------------------------------------- %% examp18.2-1 %+++++++++++++++++++++++++调用normfit函数求解+++++++++++++++++++++++++++++++ x = [15.14 14.81 15.11 15.26 15.08 15.17 15.12 14.95 15.05 14.87]; [muhat,sigmahat,muci,sigmaci] = normfit(x,0.1) %++++++++++++++++++++++++++++调用mle函数求解++++++++++++++++++++++++++++++++ x = [15.14 14.81 15.11 15.26 15.08 15.17 15.12 14.95 15.05 14.87]; [mu_sigma,mu_sigma_ci] = mle(x,'distribution','norm','alpha',0.1) %% examp18.2-2 x = [0.7917,0.8448,0.9802,0.8481,0.7627 0.9013,0.9037,0.7399,0.7843,0.8424 0.9842,0.7134,0.9959,0.6444,0.8362 0.7651,0.9341,0.6515,0.7956,0.8733]; PdfFun = @(x,theta) theta*x.^(theta-1).*(x>0 & x<1); [phat,pci] = mle(x(:),'pdf',PdfFun,'start',1) %% examp18.3-1 %++++++++++++++++++++++++++++++++双侧检验++++++++++++++++++++++++++++++++++ x = [97 102 105 112 99 103 102 94 100 95 105 98 102 100 103]; [h,p,muci,zval] = ztest(x,100,2,0.05) %++++++++++++++++++++++++++++++++单侧检验++++++++++++++++++++++++++++++++++ x = [97 102 105 112 99 103 102 94 100 95 105 98 102 100 103]; [h,p,muci,zval] = ztest(x,100,2,0.05,'right') %% examp18.3-2 x = [49.4 50.5 50.7 51.7 49.8 47.9 49.2 51.4 48.9]; [h,p,muci,stats] = ttest(x,50,0.05) %% examp18.3-3 x = [20.1, 20.0, 19.3, 20.6, 20.2, 19.9, 20.0, 19.9, 19.1, 19.9]; y = [18.6, 19.1, 20.0, 20.0, 20.0, 19.7, 19.9, 19.6, 20.2]; alpha = 0.05; tail = 'both'; vartype = 'equal'; [h,p,muci,stats] = ttest2(x,y,alpha,tail,vartype) %% examp18.3-4 x = [49.4 50.5 50.7 51.7 49.8 47.9 49.2 51.4 48.9]; var0 = 1.5; alpha = 0.05; tail = 'both'; [h,p,varci,stats] = vartest(x,var0,alpha,tail) %% examp18.3-5 x = [20.1, 20.0, 19.3, 20.6, 20.2, 19.9, 20.0, 19.9, 19.1, 19.9]; y = [18.6, 19.1, 20.0, 20.0, 20.0, 19.7, 19.9, 19.6, 20.2]; alpha = 0.05; tail = 'both'; [h,p,varci,stats] = vartest2(x,y,alpha,tail)